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ABSTRACT 

We calculate collision strengths and their thermally-averaged Maxwellian values for 
electron excitation and de-excitation between the fifteen lowest levels of singly-ionised 
cobalt, Co^, which give rise to emission lines in the near- and mid-infrared. Transition 
probabilities are also calculated and relative line intensities predicted for conditions 
typical of supernova ejecta. The diagnostic potential of the 10.52, 15.46 and 14.74 p,m 
transition lines is briefly discussed. 

Key words: atomic data - atomic processes - radiation mechanisms: non-thermal - 
supernovae: general - infrared: general. 


1 INTRODUCTION 

Cobalt transition lines are useful astronomical probes es¬ 
pecially in analysing the observational data of supernovae 
(SNe) where these emissions can be used to assess the nu¬ 
cleosynthesis and decay processes associated with the SN 
explosion. The spectra of supernova ejecta show prominent 
infrared forbidden lines of singly- and doubly-ionised ions 
of the iron group elements: nickel, cobalt and iron, where 


dioactive decay chain ®®Ni —>-®®Co —^®®Fe (Colgate & McKee 

1969| IKuchner et al\ 

19941 IBowers et a/|19971|Liu et a/|19971 

Churazov et al|2014 

. Hence, the lines of these elements can 


be used to examine the underlying nuclear processes. These 
lines can also be used for the determination of particle num¬ 
ber density and for analysing the thermodynamic conditions 
of these objects ( Nussbaun]er^&^^Storey||19^ jKuclmer^e£^ 


19941 IBowers et a/||1997 

Adelman et a'||2000| |del Peloso et 

aZ|2005 Bergemann et al 

2010^ 


Observations of cobalt lines in super and symbiotic no¬ 
vae, largely forbidden transitions in the infrared and optical 
parts of the spectrum, have been reported in a number of 


studies ( 

Axelrod 

1980 

Li et al 

1993 

Jennings et al 1993 

Kuchner et a/|1994| |Uinerstein|1995| | 

lowers et a/|19971 |Liu| 

et al|1997 

McKenna et al|1997 

Churazov et al|2014 

I. Defi- 


nite or tentative observations of certain forbidden transitions 
of Co in the spectra of planetary nebulae and H ll regions 


have also been reported (Baluteau et al 

1995 

Sharpee et 

al 20031 Esteban et al\ 20041 Zhang et al\ 

200^ 

’ottasch & 

Surendiranath||20051 Sharpee et a/||20071 Wang & Liu||20071 


E-mail: pjs@star.ucl.ac.uk 


Fang fc LiupOlIb Cobalt lines have also been observed in 


other astronomical objects like early-type stars (Adelman et 


aZ|2000l and cool stars (Bergemann et al|2010l as well as in 

the solar spectrum (|Kerola & Aller 

1976 

Salih et al 1985 

Pickering et al|1998 [Bergemann et a/|2010 

. Relevant obser- 


vational data related to the SN spectra in general and to the 
IR lines in particular have been gathered in the past using 
various instruments such as the Infrared Spectrograph (IRS) 
on the Spitzer Space Telescope and Himalaya Faint Object 
Spectrograph Camera (HFOSC) on the Himalayan Chan¬ 


dra Telescope (see e.g. Sahu et al (20061 and Jerkstrand et 


There are several experimental investigations (see e.g. 


Sugar & Corliss (1981 19851; Salih et al (19851; 

Crespo 

Lopez-Urrutia et al (19941; Mullman et al 

(1998a|bl; 

Picker- 


ing et al|(|1998[);|Ruffon &: Pickering|(|2013 1) in which atomic 


data, such as energy levels and radiative transition proba¬ 
bilities, have been reported. 

[Nussbaumer fc Storey] ( |1988[ ) reported calculations of 
forbidden transition probabilities between the energetically 
lowest eight levels of Co"*" using a multi-configuration atomic 
model. The calculations, which are based on employing the 
Superstructure code (Eissner et al|[l974 Nussbaumer & 


Storey 19781, have also included summarisation of simi¬ 


lar calculations for Co l and Co ill transitions. These data 
were intended to facilitate the analysis of the observations 
of SNe in general and SN1987A in particular. [Raassen et ai\ 
(19981 conducted experimentally-based calculations of os¬ 


cillator strengths and transition probabilities for a number 
of dipole-allowed and dipole-forbidden transitions of singly- 
ionised cobalt using a semi-empirical approach based on 
employing an orthogonal operator method within an inter- 
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2 P.J. Storey & C.J. Zeippen & T. Sochi 


mediate coupling scheme. The data were intended for use 
in analysing astronomical spectra from such objects as Co 
stars and late type supernovae. For similar purposes, |Quinet| 
119981 computed radiative transition probabilities for the 
dipole-forbidden transitions between the 47 metastable en¬ 
ergy levels in the 3d®, 3d^4s and 3d®4s^ configurations of 
the singly-ionised cobalt using a configuration interaction 
Hartree-Fock relativistic approach with optimised radial pa¬ 
rameters based on the available observed energy levels. 

In this paper, we report computational atomic data in 
the form of collision strengths and their thermally-averaged 
Maxwellian values for electron excitation and de-excitation 
between the 15 lowest levels of Co'*' as well as transition 
probabilities for a number of the forbidden lines of Co ll 
in the infrared section of the spectrum. The investigation 
is generally based on employing the R-matrix method and 


codej^ (Berrington et al 


1974 


1987 


Hummer et al 


1993 


Berrington et a/|1995 1 where the scattering calculations are 


conducted using a 13-configuration atomic target within a 
Breit-Pauli intermediate coupling approximation. 

In Section of the present paper, we discuss the tar¬ 
get used in our Co"^ model and give the resulting transition 
probabilities. Details of the Breit-Pauli R-matrix electron 
scattering calculations can be found in Section Results 
and some diagnostics appear in Section]^ The paper is fi¬ 
nalised in Sectionj^with general conclusions and discussion. 


2 Co+ ATOMIC STRUCTURE 


2.1 The scattering target 

A schematic diagram of the term structure of Co ll is shown 
in Figure[^ The lowest 19 terms are of even parity from the 
configurations 3d® and 3d^4s. In this work we are concerned 
with the excitation mechanisms of the forbidden transitions 
among these even-parity terms, particularly those lying in 
the near- and mid-infrared. The levels giving rise to mid- 
infrared transitions and the corresponding wavelengths are 
shown in Figure Transitions from higher terms give rise 
to near-infrared lines but these are expected to be weaker 
at the temperatures typical of supernova ejecta. 

It is expected that the odd-parity terms of the 3d^4p 
configuration will give rise to resonance series that affect the 
collision strengths for excitation of the low-lying even-parity 
levels and should therefore be included in the target for the 
electron scattering calculation. The extent of our target is 
shown by the solid line in Figure and includes 26 terms 
and 72 levels. The target states were expanded over the set 
of thirteen electron configurations listed in Table and the 
target wavefunctions were calculated with the program Au- 


TOSTRUCTURE, (Eissner et al|[T974 Nussbaumer & Storey 


1978 Badnell 1|2011 1, which uses radial wavefunctions cal 

culated in scaled Thomas-Fermi-Dirac statistical model po¬ 
tentials. The scaling parameters were determined by min¬ 
imising the sum of the energies of all the target terms, com¬ 
puted in LS'-coupling, i.e. neglecting all relativistic effects. 
The resulting scaling parameters are given in Table where 


1 See Badnell: R-matrix write-up on WWW. URL: 
amdpp.phys.strath.ac.uk/UK_RmaX/codes/serial/WRITEUP. 
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Figure 1. Schematic term energy diagram of Co II. The solid line 
shows the extent of the close-coupled target states. 
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Figure 2. The lowest levels of Co+ and the mid-infrared tran¬ 
sitions between them. The energy scale on the right is in cm“^ 
while the transition wavelengths next to the arrows are in /rm. 


a negative scaling parameter, A„;, signifies a hydrogenic cor¬ 
relation orbital with nuclear charge 27|A„;|. 

In Table[^we compare the term energies calculated with 
our target with experiment for the 26 terms of the target. 
The calculated term energies include one-body relativistic 
shifts, the mass and Darwin terms, and the spin-orbit in¬ 
teraction. This is the level of approximation that applies in 
the R-matrix code used for the electron scattering calcula¬ 
tion. In Table|^we compare the calculated energies of the 15 
lowest levels with experimental values. In this table we also 
show the results obtained when the two-body fine structure 
interactions described by Eissner et al (19741 are included. 
The overall fine-structure splittings of the tabulated terms 
are significantly improved when the two-body terms are in- 
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Atomic data for Co ll infrared lines 3 


Table 1. The target configuration basis where the core structure 
is suppressed. The bar indicates a correlation orbital. 


“Sc? 

3d'^ 4s 
3d® 4 s 2 
3d® 4p2 

3d'^4d 
3d'^5s 
3d'^5d 
3d'^6d 
3d® 4s Id 
3d® 4d2 
3d® 4 s 2 Id 


3d'^ 4p 
3d®4s4p 


Table 3. Energies of the 26 target terms in cm“^. The calcu- 
lated values include only the spin-orbit contribution to the fine- 
structure energies. Core structure is suppressed from all configu¬ 
rations. 


Config. Term 


3d® 

3d'^4s 

3d® 

3d'^4s 

3d® 

3d'^4s 


Table 2. Potential scaling parameters. The bar over the principal 
quantum number and the minus sign attached to the value of a 
scaling parameter signifies a correlation orbital. 


®F 

®F 

Id 

3p 

5p 

iG 

®G 

3p 

3p 

iG 

®H 

ip 

®D 

irr 


Term Energy 
Exp.t Calc. 


0 

3457 

9774 

10954 

12648 

17275 

18277 

21261 

23479 

24441 

24450 

26747 

26888 

27353 


0 

2961 

9570 

13639 

15349 

20043 

20935 

23588 

26652 

29699 

27245 

29406 

33570 

31649 


Is 

1.43187 

2p 

1.08106 

3d 

1.02832 


Id 

30502 

35335 

2s 

1.13734 

3p 

1.04220 

4d 

-0.34647 

3d’^4s 

®F 

40228 

46895 

3s 

1.05922 

4p 

0.98776 

5d 

-2.19279 

3d®4s2 

®D 

40560“ 

41967 

4s 

0.98959 



6d 

-2.00375 

3d® 

is 

42303 

51388 

5s 

-0.55000 





3d'^4s 

If 

43394 

50213 







3d'^4p 

5po 

45018 

42614 


eluded, with the average difference from experiment falling 
from 5.3% to 1.8%. Table serves as a key to the levels 
for use in later results for transition probabilities, collision 
strengths and effective collision strengths. 

A further measure of the quality of the target is a com¬ 
parison between weighted oscillator strengths, gf, calculated 
in the length and velocity formulations. Good agreement be¬ 
tween the two formulations is a necessary but not sufficient 
condition for ensuring the quality of the target wavefunc- 
tions. This comparison, given in Table shows an average 
difference in the absolute values of gf of 10.9% between the 
two formulations, which we consider acceptable for an open 
d-shell system. 


2.2 Transition probabilities 

Using the target wavefunctions described above, we com¬ 
puted the forbidden transition probabilities among the low- 
lying even parity terms. The calculated energies are replaced 
by experimental energies to correct the energy factors con¬ 
necting the ab initio calculated line strengths to the tran¬ 
sition probabilities. The results are given in Table with 
comparison to values obtained from previous investigations. 
We only tabulate those probabilities from a given upper level 
which exceed 1% of the total probability from that level. We 
also tabulate the magnetic dipole transition probabilities be¬ 
tween the levels of the and ®F terms calculated assuming 
the states are described by pure LS-coupling, in which case 


5D° 
5 GO 
3 GO 

3po 

3D° 


46319 

46905 

48507 

49518 

51289 


44049 

44577 

46823 

47843 

49868 


1 Experimental energies are from Sugar & Gorliss (1985 i. 
“The ®Do energy, which is not known experimentally, was 
derived from the calculated ®Di-®Do energy separation. 


Table 4. Energies in cm ^ of the 15 lowest target levels. 


Index 

Gonfig. 

Level 

Exp.i 

Calc.^ 

Calc.® 

1 

3d® 

a®F4 

0. 

0. 

0. 

2 


a®F3 

951 

988 

972 

3 


a®F2 

1597 

1676 

1645 

4 

3dl'4s 

a®F5 

3351 

2899 

2876 

5 


a®F4 

4029 

3562 

3550 

6 


a®F3 

4561 

4088 

4082 

7 


a®F2 

4950 

4476 

4474 

8 


a^Fi 

5205 

4732 

4731 

9 

3dl’4s 

b®F4 

9813 

9642 

9623 

10 


b®F3 

10709 

10531 

10523 

11 


b®F2 

11322 

11156 

11153 

12 

3d8 

aiD2 

11651 

14367 

14378 

13 


a®P2 

13261 

16022 

15970 

14 


a3Pi 

13404 

16094 

16094 

15 


a®Po 

13593 

16305 

16297 


the line strength is given by a simple formula | Nussbaumer 
fc Store>^|1988 h 


3d' 


The mid-infrared transitions between the levels of the 
® a®F and 3d^4s a®F terms are dominated by mag¬ 
netic dipole decays which change J by unity, with electric 
quadrupole transitions being orders of magnitude smaller. 
This leads to a step-wise decay through the levels within 
each term, as illustrated in Figure The lowest a®F level, 


^ Galculated with only spin-orbit interaction. 

® As 2 plus two-body fine-structure interactions. 


a^Fs, decays by a very weak magnetic dipole transition to 
the ground a®F 4 level. We discuss this transition in more 
detail at the end of this section. 

We find very close agreement with [Nussbaumer &:| 
Storey 119881 for transition probabilities within the a^F and 
a^’F terms, with differences less than 1% between the two 
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4 P.J. Storey & C.J. Zeippen & T. Sochi 


Table 5. Weighted LS oscillator strengths, gf, in the length and 
velocity formulations from the energetically lowest three terms. 



Transition 


gfv 

gfv 

3(F 


- Sd'^p 

^G° 

0.17 

0.21 



- 

3po 

1.34 

1.55 



- 

3D° 

0.75 

0.81 

3d'^4s 

sp 

- 3d74p 

5po 

10.75 

10.70 



- 

5 DO 

7.30 

6.62 



- 

5 GO 

14.23 

14.97 

3d'^4s 

Sp 

- 3d^4p 

3 GO 

7.51 

9.09 



- 

3po 

5.82 

6.52 



- 

3do 

4.29 

4.42 


calculations. Differences of up to 45% occur for transitions 
between the 6®F and two lower terms and within that term. 
We also find reasonable agreement with the work of 


Quinet (19981 for the transitions of interest here, the mag¬ 
netic dipole transition probabilities within the a^F and a®F 
terms. His results differ by no more than 20% from the 
present work for those transitions. Larger differences, of up 
to 60% are seen for transitions from the higher PF, a^D and 
a^P terms to the ground PF term, although the remaining 
transitions show much smaller differences, typically about 
10%. The calculation of Quinet (19981 uses a similar con¬ 
figuration basis to the present work but employs a Hartree- 
Fock approach which incorporates fitting of Slater param¬ 
eters to experimental energies. Consequently, that method 
yields calculated energies in much better agreement with ex¬ 
periment than our ab initio results. We do, however, make 
empirical corrections to the LS-Hamiltonian matrix to bring 
our final calculated level energies into good agreement with 
experiment, to ensure that the spin-orbit interactions be¬ 
tween levels of different terms are corrected for any errors 
due to incorrect term energy separations. We also use exper¬ 
imental energies in the calculation of transition probabilities 
from line strengths. 

Table[^also shows the forbidden transition probabilities 
calculated by [Raassen et al (19981, which show very good 
agreement with the present work, differing by less than 10% 
in all but three cases and agreeing within 1% for all the 
magnetic dipole transitions between the levels of the a^F 
and a®F terms. Despite this very close agreement for the in¬ 
frared transitions, the larger differences between our results 
and those of Quinet (19981 should probably be viewed as 
a measure of the uncertainty in the values of the transition 
probabilities for this rather complex ion. 

The a®F 5 -a®F 4 , which was briefly referred to above, re¬ 
quires a separate discussion. Since there are no magnetic 
dipole matrix elements between different terms, this transi¬ 
tion must proceed via an int eraction between 3d^ a^F 4 and 

(|1988|), 


3d'^4s a^F 4 . As discussed by Nussbaumer & Storey 


there is no direct spin-orbit interaction between a'^F 4 and 
a®F 4 levels; only very small two-body fine-structure terms. 
A combination of configuration interaction between a^F and 
b^F, plus spin-orbit interaction between b^F 4 and a®F 4 , is 
also excluded because the matrix element for the a^F - PF 
electrostatic interaction is also zero. In practice, the tran¬ 
sition can occur via configuration interaction between a^F 
and c^F, which is non-zero, and spin-orbit interaction be¬ 
tween PF 4 and a^F 4 , but the c^F is at relatively high energy 
(40228 cm“^) so the interaction is very weak. 

In their discussion of this transition, [Nussbaumer fc| 


Storey (19881 point out that since it is very weak, it might 


be significantly altered if the relativistic corrections to the 
magnetic dipole operator described by [Eissner fc Zeippen | 
(19811 were included. These corrections are included in the 
current work and we find that they change the a®F 5 -a^F 4 
transition probability by less than 1%. Among the four low¬ 
est terms, the largest changes due to these corrections, rang¬ 
ing from 0 to 3.2%, are seen in the b^F-a^F multiplet. We 
also find that the exclusion of two-body fine-structure terms 
only changes the a®F 5 -a®F 4 transition probability by a small 
amount (5%). 

For this transition, our calculated probability of 
9.93x10“® s“^ is a factor 2.0 times larger than the value of 


Nussbaumer & Storey (19881 and a factor 3.2 times larger 


than that found by Quinet (19981. Raassen et al (19981 do 
not give a value for this transition probability. In addition 
to the mechanism described above, we find that the a®F 4 
and a®F 4 levels can interact via ®F terms belonging to the 
3d^4d configuration. We believe that it is presence of these 
interactions that causes the large differences between the 
different calculations. The present calculation and that of 
[Nussbaumer fc Storey] ( [1988P utilise a short-range correla¬ 
tion 4d orbital to approximate the effect of the infinite series 
of bound and continuum d orbitals that should be accounted 
for. [Quinet] ( |199^ does not use such orbitals and so underes¬ 
timates the contribution of the 3d^nd configurations to the 
transition probability. 

Within the a®F and a®F terms, radiative decays occur 
stepwise between adjacent levels with probabilities of order 
10“^ s“^, while the probability of a transition between the 
a®F and a®F terms is typically four orders of magnitude 
smaller. Hence the only transition of significance between 
the two terms is a®F 5 -a®F 4 . In a physical situation where 
levels are excited by electron collisions and decay by radia¬ 
tive decay or collisional de-excitation this transition alone 
determines the critical electron density above which the rate 
of collisional de-excitation of the a®F levels is larger than ra¬ 
diative decay. We return to this point in Section[4.1| 


3 SCATTERING CALCULATIONS 


The Breit-Pauli R-matrix method which is used in this cal¬ 


culation is described fully elsewhere (Hummer et al 1993 
Herrington et a(|[1995 1 and references therein. The calcula¬ 


tions described here were made with the parallel versions of 
the codes|^ We use an R-matrix boundary radius of 14.9 au, 
to encompass the most extended target orbital (4p). The 
expansion of each scattered electron partial wave is over a 
basis of 12 functions within the R-matrix boundary, and the 
partial wave expansion extends to a maximum of 2J = 19. 

Collision strengths are computed at 20000 equally 
spaced values of the energy in the resonance region and a fur¬ 
ther 1000 values in the region where all scattering channels 
are open, up to an incident electron energy of 1.5 Rydberg. 
In Figurej^we illustrate our results with the calculated colli¬ 
sion strengths between the three levels of the ground 3d® ®F 
term as a function of final electron energy up to 0.4 Rydberg 


^ See Badnell: R-matrix write-up on WWW. URL: 
http://amdpp. phys.strath.ac.uk/tamoc/codes/parallel/PWRITEUP. 
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O 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



O 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Figure 3. Collision strength (vertical axis) versus final electron 
energy in Rydberg (horizontal axis) for the (a) 1-2, (b) 1-3 and 
(c) 2-3 transitions. For level indexing refer to Table 


above threshold. Due to the multiple close lying thresholds 
the collision strengths display very dense and complex reso¬ 
nance structures. 


4 RESULTS AND DISCUSSION 

In Table the final thermally averaged collision strengths 
between the 15 energetically lowest levels are given as a func¬ 
tion of electron temperature. In the energy region where all 
scattering channels are open we find some small irregular 
features in the collision strengths that are almost certainly 
non-physical and caused by the correlation orbitals in the 
target representation. We have computed thermally aver¬ 
aged collision strengths for the transitions and temperature 
range given in Tablej^both including and excluding the con¬ 
tribution from the region of all channels open, and find the 
largest change for any transition is 0.3% at logioT = 4.0, 
2.4% at logioT = 4.2 and 9.4% at logioT = 4.4. For the 
excitations from the a^F levels to the a®F levels, the max¬ 
imum difference at any temperature is 2.2%. Hence for the 
infrared transitions of interest here the contribution from 
any non-physical features is insignificant. The values tabu¬ 
lated in Table were computed using the full energy range. 


4.1 The IR line ratios 

The transitions between the levels of the a^F and a®F 
terms give rise to mid-infrared lines lying between 10.5 and 
39.3 ^m. The relative intensity of the two lines arising from 



Figure 4. Emissivity ratios of pairs of Co II infrared lines, show¬ 
ing contours of constant electron temperature (dashed lines) and 
electron number density (solid lines). Lines are labelled with the 
logio of temperature in K or density in cm“^. 


the a®F term (10.52 ^m and 15.46 ^m, see Figure is 
sensitive to electron density for densities less than about 
10® cm“® and relatively insensitive to temperature for tem¬ 
peratures greater than about 3000 K due to the low exci¬ 
tation energies. The lines from within the higher lying a®F 
term, the 14.74 fim line for example, are more sensitive to 
electron temperature. Ratios of the intensities of these three 
lines can therefore be used to simultaneously determine the 
electron temperature and density of the emitting material if 
they lie in the appropriate range. This is illustrated in Fig¬ 
ure where we plot the 10.52/14.74 /rm ratio against the 
10.52/15.46 fj,m ratio for va rious temperatures and densi¬ 
ties. We discussed in Section [2.2| the role of the a®F 5 -a®F 4 
transition in determining the relative populations of the a®F 
and a®F levels and the rather large differences in the cal¬ 
culated values for the radiative decay probability for this 


transition. The value calculated by Quinet (19981 is about 


a factor of three smaller than ours and a lower transition 
probability would mean that the critical density for colli- 
sional de-excitation of the a®F levels would be three times 
smaller. Above the critical density the relative level popu¬ 
lations tend to their Boltzmann values and depend only on 
temperature as shown by the lines in Figure 1^ correspond¬ 
ing to electron number densities of 10^ and 10^ cm“®. If the 


transition probability of Quinet (19981 were adopted, these 


high density limit values of the ratios would be reached at 
a density three times lower. 


5 CONCLUSIONS 

In this paper, the Co ll forbidden transitions between the fif¬ 
teen lowest energy levels of singly-ionised cobalt, Co^, have 
been investigated. Radiative transition probabilities and, for 
the first time for this ion, collision strengths for excitation 
by electrons and their thermally-averaged values based on a 
Maxwell-Boltzmann electron energy distribution have been 
calculated. An elaborate atomic target was used to perform 
the scattering calculations using the R-matrix method in 
the Breit-Pauli approximation and intermediate coupling 
scheme. 
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6 P.J. Storey & C.J. Zeippen & T. Sochi 


We have compared the radiative transition probabili¬ 
ties with those from three previons calculations by |Nuss- 


baumer & Storey (19881, Quinet (19981 and Raassen et al 


(19981, and briefly discussed the consequences of the dif¬ 
ferences seen for the spectroscopy of the infrared lines. Al¬ 
though there is good agreement between our data and the 
previous works, there are also some significant differences, 
particularly for the important a®F 5 -a®F 4 transition. This 
probability is clearly very sensitive to the choice of configu¬ 
rations and orbitals, and we have argued that our basis and 
orbitals give the most reliable value for this probability. 

The diagnostic potential of the infrared lines is illus¬ 
trated by showing that the three lines at 10.52, 14.74 and 
15.46 pm can potentially be used to measure the electron 
density and temperature of the emitting region. These lines 
are of particular use in the study of supernovae. 


6 ACKNOWLEDGMENTS 

PJS and CJZ acknowledge financial support from the 
Atomic Physics for Astrophysics Project (APAP) funded 
by the Science and Technology Facilities Council (STFC). 
The present collaboration benefited from visits to Meudon 
by PJS from 1998 to 2013, with support provided by the 
Observatoire de Paris and by the Universite Paris 7. The 
hospitality of the Observatoire de Paris was much appreci¬ 
ated. 


REFERENCES 

Adelman S.J., Gulliver A.F., Loden, L.O., 2000, A&A, 353, 
335 

Axelrod T.S., 1980, PhD thesis (UCRL-52994), University 
of California Santa Cruz 

Badnell, N.R., 2011, Comput. Phys. Commun. 182, 1528 

Baluteau J.P., Zavagno A., Morisset C., Pequignot D., 
1995, A&A, 303, 175 

Bergemann M., Pickering J.C., Gehren T., 2010, MNRAS 
401, 1334 

Berrington K.A., Burke P.G., Chang J.J., Chivers A.T., 
Robb W.D., Taylor K.T., 1974, Comp. Phys. Comm., 8, 
149 

Berrington K.A., Burke P.G., Butler K., Seaton M.J., 
Storey P.J., Taylor K.T., Yu Yan., 1987, J. Phys. B, 20, 
6379 

Berrington K.A., Eissner W.B., Norrington P.H., 1995, 
Comp. Phys. Comm., 92, 290 

Bowers E.J.C., Meikle W.P.S., Geballe T.R. et al, 1997, 
MNRAS, 290, 663 

Churazov E., Sunyaev R., Isern J. et al, 2014, Nature, 512, 
406 

Colgate S.A., McKee C., 1969, ApJ, 157, 623 

Crespo Lopez-Urrutia J.R., Ulbe M., Neger T., Jager H., 
1994, J. Quant. Spect. Rad. Trans., 52, 89 

Dinerstein H.L., 1995, “Infrared emission lines as probes of 
gaseous nebulae”, in: “The Analysis of Emission Lines”, 
eds. Williams R.E. and Livio M., Space Telescope Science 
Institute Symposium Series 8 

Eissner W., Jones M. & Nussbaumer H., 1974, Comput. 
Phys. Commun. 8, 270 


Eissner, W., Zeippen, C.J., 1981, J. Phys. B. 14, 2125 
Esteban C., Peimbert M., Garcfa-Rojas J., Ruiz M.T., Pe- 
imbert A., Rodriguez M., 2004, MNRAS, 355, 229 
Fang X., Liu X.-W., 2011, MNRAS, 415, 181 
Hummer D.G., Berrington K.A., Eissner W., Pradhan 
A.K., Saraph H.E., Tully J.A., 1993, A&A, 279, 298 
Jennings D.E., Boyle R.J., Wiedemann G.R., Moseley S.H., 

1993, ApJ, 408, 277 

Jerkstrand A., Fransson C., Maguire K., Smart! S., Ergon 
M., Spyromilio J., 2012, A&A, 546, A28 
Kerola D.X., Aller L.H., 1976, Publications of the Astro¬ 
nomical Society of the Pacific, 88, 122 
Kuchner M.J., Kirshner R.P., Pinto P.A., Leibundgut B., 

1994, ApJ, 426, L89 

Li H., McCray R., Sunyaev R.A., 1993, ApJ, 419, 824 
Liu W., Jeffery D.J., Schultz D.R., Quinet P., Shaw J., 
Pindzola M.S., 1997, ApJ, 489, L141 
McKenna F.C., Keenan F.P., Hambly N.C. et al, 1997, 
ApJS, 109, 225 

Mullman K.L., Cooper J.C., Lawler J.E., 1998, ApJ, 495, 
503 

Mullman K.L., Lawler J.E., Zsargo J., Federman S.R., 
1998, ApJ, 500, 1064 

Nussbaumer H., Storey P.J., 1978, A&A, 64, 139 
Nussbaumer H, Storey P.J., 1988, A&A, 200, L25 
del Peloso E.F., Cunha K., da Silva L., Porto de Mello G.F., 
2005, A&A, 441, 1149 

Pickering J.C., Raassen A.J.J., Uylings P.H.M., Johansson 
S., 1998, ApJS, 117, 261 

Pottasch S.R., Surendiranath R., 2005, A&A, 444, 861 
Quinet P., 1998, A&AS, 129, 147 

Raassen A.J.J., Pickering J.C., Uylings P.H.M., 1998, 
A&AS 130, 541 

Ruffon M.P., Pickering J.C., 2013, ApJS, 207, 20 
Sahu D.K., Anupama G.C., Srividya S., Muneer S., 2006, 
MNRAS, 372, 1315 

Salih S., Lawler J.E., Whaling W., 1985, Phys. Rev. A, 31, 
744 

Sharpee B., Williams R., Baldwin J.A., van Hoof P.A.M., 
2003, ApJS, 149, 157 

Sharpee B., Zhang Y., Williams R., Pellegrini E., Cav- 
agnolo K., Baldwin J.A., Phillips M., X.-W. Liu, 2007, 
ApJ, 659, 1265 

Sugar J., Corliss C., 1981, J. Phys. Chem. Ref. Data, 10, 
1097 

Sugar J., Corliss C., 1985, J. Phys. Chem. Ref. Data, 14, 1 
Wang W., Liu X.-W., 2007, MNRAS, 381, 669 
Zhang Y., Liu X.-W., Luo S.-G., Pequignot D., Barlow 
M.J., 2005, A&A, 442, 249 


© 0000 RAS, MNRAS 000, 000-000 

















Atomic data for Co ll infrared lines 7 


Table 6. Transition probabilities in s“^ among the energetically lowest 15 leve ls (see Table for level indexing) as obtained from the 
current work (CW), [Nussbaumer & Storey] i flQSSl l (NS), [Quinetj ( [l998[ | (Q) and |Raassen et aZ| l |199^ (RPU) where i and j refer to the 
lower and upper levels respectively. Only transition probabilities that are at least 1% of the total probability from a given upper level 
are listed. The powers of 10 by which the number is to be multiplied are given in brackets. 


Transition 

j ^ 

CW 

NS 

Q 

RPU 


Transition 

j 

CW 

i 

NS 

Q 

RPU 

2 

1 

2.23(-2) 

2.23(-2) 

2.21(-2) 

2.24(-2) 

10 

9 

1.88(-2) 

1.87(-2) 

1.67(-2) 

1.88(-2) 

3 

2 

9.71(-3) 

9.73(-3) 

9.88^3) 

9.73^3) 

11 

2 

2.73^2) 

1.87^2) 

3.16(-2) 

1.72^2) 

4 

1 

9.13(-6) 

4.90(-6) 

3.09^6) 

— 

11 

3 

2.81^2) 

3.09(-2) 

4.06(-2) 

2.57(-2) 

5 

4 

1.24^2) 

1.23(-2) 

1.03^2) 

1.24(-2) 

11 

7 

4.61^3) 

5.30^3) 

4.34(-3) 

4.39^3) 

6 

5 

1.08^2) 

1.08(-2) 

9.27(-3) 

1.09^2) 

11 

8 

1.41(-2) 

1.69(-2) 

1.34(-2) 

1.35^2) 

7 

6 

5.72(-3) 

5.72 (-3) 

5.00^3) 

5.74(-3) 

11 

10 

7.82(-3) 

8.32(-3) 

7.26(-3) 

8.19^3) 

8 

3 

4.33(-5) 

— 

— 

— 

12 

1 

7.61^3) 

— 

9.74(-3) 

— 

8 

7 

1.78(-3) 

1.78(-3) 

1.57(-3) 

1.79(-3) 

12 

2 

1.67(-1) 

— 

1.74(-1) 

— 

9 

1 

3.36^2) 

3 . 85 I 2 ) 

5.36^2) 

3.43^2) 

12 

3 

7.99^2) 

— 

8.16^2) 

1.27(-3) 

9 

2 

3.98^3) 

4.91(-3) 

6.33^3) 

4.01^3) 

13 

1 

3.34(-2) 

— 

4.98(-2) 

— 

9 

4 

2.90^2) 

2 . 89 I 2 ) 

2.66^2) 

2.81^2) 

13 

2 

6.98^2) 

— 

6.70^2) 

— 

9 

5 

3.27^3) 

2.82(-3) 

2.98(-3) 

3.07^3) 

13 

3 

1.82(-2) 

— 

1.63(-2) 

— 

9 

6 

5.21^3) 

4 . 28 I 3 ) 

4.81^3) 

5.05^3) 

13 

12 

3.01^2) 

— 

2.96(-2) 

— 

10 

1 

1.45^2) 

1.63(-3) 

2.30^2) 

1.51^2) 

14 

2 

2.92(-2) 

— 

4.20(-2) 

— 

10 

2 

2.22^2) 

2 . 58 I 2 ) 

3.47^2) 

2.25^2) 

14 

3 

1.19^2) 

— 

1.68^2) 

— 

10 

3 

5.60^3) 

6.97(-3) 

8.74(-3) 

5.65(-3) 

14 

11 

1.99^3) 

— 

— 

— 

10 

5 

4.43^3) 

6.41(-3) 

4.17^3) 

3.97^3) 

14 

12 

2.78 (-2) 

— 

2.55(-2) 

— 

10 

6 

4.95^3) 

4.91(-3) 

4.53^3) 

4.57^3) 

15 

3 

3.70(-2) 

— 

5.34(-2) 

— 

10 

7 

1.04(-2) 

I.O 4 I 2 ) 

9.59^3) 

9.65^3) 
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Table 7: Thermally-averaged collision strengths among the 15 energetically lowest levels as a function of logic of temperature 
in K where i and j refer to the index of the lower and upper level respectively (see Table ^ for indexing). 


* j logioT 




2.0 

2.2 

2.4 

2.6 

2.8 

3.0 

3.2 

3.4 

3.6 

3.8 

4.0 

4.2 

4.4 

1 

2 

4.250 

4.079 

3.969 

3.899 

3.760 

3.519 

3.308 

3.297 

3.412 

3.420 

3.235 

2.925 

2.563 

1 

3 

1.449 

1.320 

1.189 

1.064 

0.947 

0.845 

0.781 

0.775 

0.802 

0.807 

0.772 

0.709 

0.636 

1 

4 

3.202 

3.326 

3.404 

3.400 

3.369 

3.405 

3.512 

3.579 

3.495 

3.242 

2.891 

2.537 

2.209 

1 

5 

2.292 

2.136 

1.988 

1.918 

1.978 

2.138 

2.285 

2.314 

2.191 

1.955 

1.684 

1.447 

1.251 

1 

6 

0.674 

0.751 

0.834 

0.937 

1.054 

1.149 

1.186 

1.150 

1.049 

0.909 

0.767 

0.656 

0.572 

1 

7 

0.235 

0.241 

0.271 

0.319 

0.366 

0.395 

0.401 

0.385 

0.349 

0.301 

0.254 

0.220 

0.196 

1 

8 

0.037 

0.039 

0.046 

0.057 

0.071 

0.081 

0.084 

0.081 

0.073 

0.063 

0.053 

0.047 

0.043 

1 

9 

2.091 

2.135 

2.169 

2.196 

2.212 

2.213 

2.192 

2.142 

2.065 

1.981 

1.941 

1.975 

2.004 

1 

10 

0.596 

0.589 

0.586 

0.591 

0.607 

0.625 

0.633 

0.626 

0.609 

0.591 

0.594 

0.630 

0.662 

1 

11 

0.143 

0.146 

0.150 

0.156 

0.161 

0.161 

0.155 

0.146 

0.136 

0.127 

0.128 

0.139 

0.149 

1 

12 

0.679 

0.726 

0.809 

0.926 

1.033 

1.084 

1.080 

1.055 

1.029 

1.003 

0.973 

0.942 

0.915 

1 

13 

0.641 

0.654 

0.654 

0.637 

0.618 

0.620 

0.642 

0.672 

0.706 

0.741 

0.769 

0.783 

0.773 

1 

14 

0.253 

0.251 

0.241 

0.228 

0.218 

0.216 

0.220 

0.228 

0.243 

0.268 

0.297 

0.320 

0.328 

1 

15 

0.074 

0.070 

0.065 

0.062 

0.061 

0.063 

0.066 

0.069 

0.072 

0.077 

0.086 

0.093 

0.095 

2 

3 

4.806 

4.372 

3.886 

3.409 

2.967 

2.590 

2.361 

2.357 

2.475 

2.513 

2.399 

2.182 

1.919 

2 

4 

1.081 

1.087 

1.068 

1.026 

0.977 

0.946 

0.936 

0.924 

0.880 

0.799 

0.700 

0.608 

0.528 

2 

5 

2.013 

2.021 

1.980 

1.902 

1.828 

1.790 

1.783 

1.762 

1.680 

1.527 

1.337 

1.157 

0.997 

2 

6 

1.445 

1.437 

1.434 

1.435 

1.456 

1.500 

1.543 

1.541 

1.463 

1.318 

1.147 

0.993 

0.863 

2 

7 

0.912 

0.868 

0.837 

0.839 

0.875 

0.926 

0.963 

0.960 

0.906 

0.811 

0.704 

0.612 

0.538 

2 

8 

0.258 

0.270 

0.296 

0.340 

0.390 

0.431 

0.453 

0.452 

0.426 

0.381 

0.331 

0.290 

0.257 

2 

9 

0.733 

0.795 

0.835 

0.847 

0.836 

0.810 

0.778 

0.743 

0.708 

0.674 

0.654 

0.654 

0.649 

2 

10 

1.314 

1.334 

1.317 

1.273 

1.225 

1.182 

1.137 

1.086 

1.031 

0.979 

0.957 

0.980 

1.006 

2 

11 

0.486 

0.484 

0.490 

0.509 

0.531 

0.545 

0.544 

0.531 

0.512 

0.495 

0.495 

0.518 

0.540 

2 

12 

0.448 

0.455 

0.497 

0.579 

0.672 

0.732 

0.749 

0.740 

0.723 

0.700 

0.676 

0.649 

0.617 

2 

13 

0.550 

0.554 

0.551 

0.544 

0.538 

0.544 

0.562 

0.581 

0.596 

0.607 

0.616 

0.619 

0.610 

2 

14 

0.370 

0.397 

0.409 

0.406 

0.397 

0.391 

0.389 

0.388 

0.392 

0.402 

0.418 

0.433 

0.436 

2 

15 

0.070 

0.064 

0.058 

0.054 

0.053 

0.055 

0.059 

0.065 

0.073 

0.084 

0.094 

0.100 

0.103 

3 

4 

0.173 

0.166 

0.158 

0.149 

0.141 

0.137 

0.139 

0.139 

0.132 

0.118 

0.101 

0.086 

0.075 

3 

5 

0.783 

0.771 

0.747 

0.712 

0.673 

0.642 

0.623 

0.606 

0.574 

0.520 

0.455 

0.394 

0.341 

3 

6 

1.197 

1.183 

1.151 

1.100 

1.051 

1.028 

1.030 

1.030 

0.992 

0.908 

0.801 

0.698 

0.606 

3 

7 

1.208 

1.174 

1.121 

1.072 

1.054 

1.077 

1.117 

1.131 

1.090 

0.995 

0.876 

0.765 

0.669 

3 

8 

0.575 

0.577 

0.592 

0.629 

0.693 

0.769 

0.829 

0.848 

0.815 

0.742 

0.652 

0.571 

0.501 

3 

9 

0.334 

0.342 

0.338 

0.321 

0.294 

0.263 

0.233 

0.207 

0.185 

0.167 

0.156 

0.152 

0.150 

3 

10 

0.748 

0.770 

0.768 

0.743 

0.707 

0.672 

0.638 

0.605 

0.572 

0.544 

0.530 

0.533 

0.533 

3 

11 

0.839 

0.827 

0.824 

0.840 

0.869 

0.891 

0.892 

0.876 

0.850 

0.820 

0.811 

0.838 

0.864 

3 

12 

0.265 

0.283 

0.318 

0.381 

0.453 

0.501 

0.515 

0.510 

0.497 

0.483 

0.475 

0.465 

0.441 

3 

13 

0.465 

0.467 

0.453 

0.433 

0.417 

0.413 

0.422 

0.434 

0.442 

0.448 

0.454 

0.457 

0.452 

3 

14 

0.254 

0.262 

0.272 

0.279 

0.283 

0.284 

0.285 

0.284 

0.288 

0.298 

0.311 

0.320 

0.321 

3 

15 

0.178 

0.168 

0.158 

0.149 

0.143 

0.141 

0.141 

0.141 

0.143 

0.147 

0.152 

0.158 

0.160 

4 

5 

9.862 

9.916 

9.923 

9.890 

9.662 

9.214 

8.727 

8.312 

7.876 

7.307 

6.712 

6.308 

5.996 

4 

6 

0.395 

0.419 

0.485 

0.581 

0.666 

0.703 

0.686 

0.630 

0.549 

0.469 

0.440 

0.499 

0.594 

4 

7 

0.127 

0.156 

0.196 

0.239 

0.274 

0.285 

0.271 

0.240 

0.202 

0.167 

0.150 

0.162 

0.187 

4 

8 

0.073 

0.075 

0.082 

0.095 

0.105 

0.104 

0.093 

0.077 

0.061 

0.049 

0.042 

0.045 

0.051 

4 

9 

18.971 

18.883 

18.780 

18.640 

18.447 

18.175 

17.721 

16.898 

15.540 

13.714 

11.799 

10.115 

8.614 

4 

10 

0.276 

0.240 

0.211 

0.190 

0.174 

0.161 

0.150 

0.145 

0.145 

0.157 

0.205 

0.294 

0.375 

4 

11 

0.044 

0.052 

0.056 

0.053 

0.047 

0.041 

0.036 

0.034 

0.034 

0.036 

0.046 

0.066 

0.085 

4 

12 

0.298 

0.297 

0.301 

0.304 

0.303 

0.297 

0.286 

0.272 

0.258 

0.244 

0.230 

0.218 

0.201 

4 

13 

0.556 

0.548 

0.549 

0.554 

0.558 

0.558 

0.553 

0.545 

0.534 

0.521 

0.508 

0.494 

0.468 

4 

14 

0.017 

0.021 

0.025 

0.027 

0.027 

0.027 

0.025 

0.024 

0.023 

0.023 

0.027 

0.033 

0.037 

4 

15 

0.036 

0.029 

0.023 

0.018 

0.014 

0.011 

0.010 

0.009 

0.008 

0.007 

0.008 

0.009 

0.009 

5 

6 

7.160 

7.608 

8.160 

8.611 

8.810 

8.785 

8.694 

8.618 

8.446 

8.035 

7.446 

6.867 

6.285 

5 

7 

0.499 

0.491 

0.506 

0.548 

0.587 

0.596 

0.570 

0.521 

0.456 

0.391 

0.368 

0.419 

0.499 

5 

8 

0.138 

0.168 

0.206 

0.240 

0.256 

0.251 

0.232 

0.205 

0.175 

0.148 

0.136 

0.151 

0.178 

5 

9 

8.909 

8.860 

8.811 

8.728 

8.556 

8.267 

7.865 

7.337 

6.647 

5.831 

5.067 

4.493 

4.005 

5 

10 

8.662 

8.640 

8.655 

8.684 

8.690 

8.634 

8.451 

8.055 

7.388 

6.496 

5.563 

4.743 

4.022 

5 

11 

0.108 

0.110 

0.110 

0.106 

0.099 

0.090 

0.083 

0.080 

0.080 

0.087 

0.117 

0.173 

0.226 

5 

12 

0.262 

0.251 

0.243 

0.234 

0.226 

0.216 

0.207 

0.197 

0.187 

0.177 

0.168 

0.161 

0.150 

5 

13 

0.357 

0.346 

0.338 

0.332 

0.327 

0.323 

0.318 

0.312 

0.303 

0.296 

0.296 

0.301 

0.296 

5 

14 

0.275 

0.276 

0.277 

0.275 

0.273 

0.270 

0.267 

0.262 

0.255 

0.248 

0.242 

0.239 

0.229 

5 

15 

0.020 

0.016 

0.013 

0.010 

0.009 

0.009 

0.009 

0.009 

0.008 

0.009 

0.010 

0.012 

0.013 

6 

7 

5.975 

6.146 

6.457 

6.876 

7.209 

7.337 

7.361 

7.372 

7.285 

6.978 

6.489 

5.970 

5.433 

6 

8 

0.360 

0.382 

0.423 

0.466 

0.489 

0.481 

0.449 

0.403 

0.350 

0.300 

0.283 

0.323 

0.386 
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Table 7: continued. 


* j logioT 




2.0 

2.2 

2.4 

2.6 

2.8 

3.0 

3.2 

3.4 

3.6 

3.8 

4.0 

4.2 

4.4 

6 

9 

2.461 

2.482 

2.488 

2.472 

2.411 

2.293 

2.131 

1.939 

1.723 

1.503 

1.342 

1.280 

1.243 

6 

10 

9.081 

9.062 

9.032 

8.955 

8.811 

8.598 

8.290 

7.826 

7.143 

6.270 

5.377 

4.611 

3.940 

6 

11 

3.168 

3.171 

3.167 

3.161 

3.148 

3.111 

3.027 

2.872 

2.629 

2.318 

2.011 

1.765 

1.552 

6 

12 

0.282 

0.254 

0.232 

0.211 

0.191 

0.174 

0.161 

0.150 

0.140 

0.129 

0.121 

0.113 

0.104 

6 

13 

0.232 

0.235 

0.232 

0.226 

0.218 

0.211 

0.203 

0.193 

0.182 

0.172 

0.171 

0.174 

0.172 

6 

14 

0.309 

0.323 

0.331 

0.329 

0.320 

0.310 

0.301 

0.291 

0.279 

0.268 

0.260 

0.256 

0.247 

6 

15 

0.076 

0.078 

0.076 

0.073 

0.071 

0.070 

0.069 

0.068 

0.066 
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